Functional renormalization flow and dynamical chiral symmetry breaking of QCD 
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The dependence of function renormalization group equation on regulators is investigated. A 
parameter is introduced to control the suppression of regulators. Functional renormalization group 
equations will become regulator-independent if this newly introduced parameter is sent to infinity 
in the end of calculation. One-loop renormalization flow equations of QCD are derived. The novelty 
is that both the coupling running equation and the mass running equation are mass-dependent. 
Different flow patterns are explored. A mechanism for non-occurrence of dynamical chiral symmetry 
breaking is arrived at. The existence of a conformal window is also discussed in the language of 
renormalization flow. 

PACS numbers: 



I. INTRODUCTION 

Renormalization is a concept about connecting physics 
at different scales. Deriving large scale physics from small 
scale physics is an important subject, and vice versa. To 
derive large scale behaviors of QCD from first principle 
is a long-standing difficult problem. 

Functional renormalization group method originates 
from Wilson's idea of integrating momentum slices suc- 
cessively, and studies how theories flow between different 
scales. The most commonly used formulations of func- 
tional renormalization group include Polchinski's formal- 
ism [l| and Wetterich's formalism [2| (for reviews, see 
for example, [!, IH). Functional renormalization group 
method uses a regulator term to modify propagators, 
thereby controls contributions of modes of different en- 
ergy scales to the effective action. In the Wetterich's 
formalism, for example, the regulator suppresses the low 
energy modes and lets the high energy modes intact. 

There have been many applications of functional re- 
mormalization group method, for example, in the anal- 
ysis of QED [5j , Non- Abeliangauge theory [g, |jj , phase 
transitions in gauge theories [8] , non-linear diffusion sys- 



tem Dyson's hierarchical model 
of applications in specific areas, see 
corssover, [HI for quark confinement 
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force, [3| for general. Recent applications include, 
sine-Gordon type models [HI, [l6| , nonrelativistic inverse 
square potential [l7j . anharmonic oscillator 18], Ising 
model [19j . etc. 

However, the use of the functional renormalization 
group method is handicapped by one issue. The flow 
equations derived through this method depend on regu- 
lation schemes, thus their physical meanings are ambigu- 
ous [20U22I ]. One bypasses this problem by arguing that 
flow equations become the true physical flow equations 
when the cutting scale runs to zero (in the Wetterich's 
formalism). But one may want not only the behavior 
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at zero energy scale, but also properties at finite energy 
scales. 

In this article, we will investigate this problem closely 
to see whether we can get a regulator-independent renor- 
malization flow. The idea is to introduce a parameter 
to control the suppressiveness of a regulator. To make 
a complete cut, the suppressiveness below the cut-off 
should be sent to infinity in the end. It is conceivable 
that the result would be regulator-independent in the in- 
finite suppression limit, since all proper regulators then 
become the same: a clear cut. 

The model that we will use is QCD in tf-dimensional 
Euclidian spacetime. Recently, lattice simulations, see 
[23M27I ] for example, show that there exists a confor- 
mal window for appropriate number of fcrmions of QCD. 
Contrary to former anticipations, dynamical chiral sym- 
metry breaking ceases to occur with a lower fermion num- 
ber than the critical fermion number for asymptotic free- 
dom, Nj — 11 x 3/2. The critical fermion number for 
chiral symmetry breaking is accessed to be in the inter- 
val 9 < Nj < 13, which is slightly larger than the lower 
turning point obtained from the two-loop /3-function of 
QCD, 51 x 3/19 « 8. 

There have been many tries to explain the existence 
of this conformal window, for example, by an ansatz of 
all- loop /3-function |28t. mass-dependent /3-functions [2!| , 
critical scaling laws 30], etc. 

In this article, we will tackle the problem through 
renormalization flows in the fh 2 -g 2 plane. Here, rh 2 
and g 2 stand for dimensionless mass and coupling, re- 
spectively. We will see that dynamical chiral symmetry 
breaking can be ceased by and only by an IR-attractive 
nontrivial fixed point with a finite dimensionless mass 
in the rh 2 -g 2 plane. All flows attracted to such a fixed 
point will have a constant dimensionless mass in the in- 
frared, which means a zero mass in the infrared. Such 
a fixed point can also exist when there is a nontrivial 
fixed point in the axis of the coupling of QCD. So, when 
Nf < 11x3/2, dynamical chiral symmetry breaking may 
not occur and a conformal window can open up. 

The structure of this articles is as follows. In section 
IIIl we briefly review the functional renormalization group 
method in Wetterich formalism; in section UITl we intro- 
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duce a parameter to control the suppressiveness of a regu- 
lator; in section lTVl we apply the method to QCD and de- 
rive its mass-dependent one-loop renormalization equa- 
tions; in section [Vj we prove the regulator-independence 
of flow equations when the suppressiveness tends to infin- 
ity; in section WT\ we draw flow patterns to gain insights 
for dynamical chiral symmetry breaking; in section I VIII 
we show that a conformal window can indeed turn up; 
finally, in section [Villi we give our conclusion. 

II. BRIEF REVIEW OF THE FUNCTIONAL 
RENORMALIZATION GROUP 

Now we sketch the functional renormalization group 
method in the Wetterich formalism. To start with, one 
adds a regulator term ASk to any given action S. For a 
scaler theory, the regulator will have the form, 

AS k = ^Jy{-q)R k {q)<p{q), (1) 

where Rk(q) = r(q 2 /k 2 )q 2 . It is to suppress the low en- 
ergy modes in the momentum integration of the effective 
action, and keep the high energy modes intact. Perform- 
ing the integration, one is left with a theory of low en- 
ergy modes. In other words, a theory of all range energy 
modes is reduced to a low energy effective theory. 

For this purpose, there are some requirements on the 
regulator. Let k be the scale at which one wants to cut off 
the theory. One first requires that Rk(q) ~ when q 2 ^> 
k 2 . Secondly, one usually takes Rk(q) ~ k 2 for q 2 <C k 2 . 
Effectively, the low energy modes gets a mass suppression 
of mass k 2 . (Later we will see that this requirement is not 
enough, which is the source of the regulator-dependence 
of functional renormalization flow.) To ensure the new 
theory identical to the original one when the regulator 
term is removed, one needs also to have Rk(q) ~ when 
k 2 ^0. 

Shown in FIG. [I] is the profile of a typical regulator. 
Plotted also is d t Rk (dt = kdk), which has a peak at 
q 2 /k 2 = 1. The peak reflects the fact that the theory is 
cut at q 2 = k 2 . 

Regulators often used may include exponentials, such 

as 

r(g2 /^ = l-e-^ e " g2/fc2 ' W 

or step-functions, such as, 

r(q 2 /k 2 ) = (k 2 /q 2 - 1)6(1 - q 2 /k 2 ), (3) 

where <d(x) is the Heaviside step function. 

Let Tk be the /c-dependent effective action of the low 
energy effective theory and T be the effective action of the 
regulated theory of action S + ASk ■ The flow equation 
of the effective action (the Wetterich equation) can be 
derived as, 

d t T k = ^Tr[(d t R k )G w ], (4) 
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FIG. 1: A typical profile of regulators. 

where <9 t = kd k - G w = (ipip) - (ip) (ip) = (r£ 2) +i? fc ) _1 , is 
the connected two-point Green function of the regulated 
theory. 

In this equation, one observes that dtTk depends 
mainly on modes where dtRk peaks. It reflects the fact 
that the renormalization flow is driven locally and de- 
pends only on modes in the vicinity of the cutting scale. 
Actually this implements Wilson's idea of 'integrating a 
single momentum slice'. 

This equation is a functional equation of Tk . It is exact 
and incorporates all non-perturbative effects. However, 
it usually cannot be solved exactly. To extract results 
from this equation, approximations must be performed. 
A common practice is the method of truncation. In this 
approach, an ansatz of the form of the effective action 
Tk is made and the right hand side of the equation (£[]) 
is projected onto the ansatz of the effective action. Then 
this equation is used to derive the flowing of coefficients 
in the ansatz of the effective action. 



III. REGULATORS 

In this section, a parameter will be introduced to con- 
trol the suppressiveness of regulators. 

As has been mentioned above, one usually requires 
Rk (q) ~ k 2 for q 2 <C k 2 . The suppression with a mass k 2 
may not be enough. Instead, we would like to replace it 
with Rk(q) ~ s 2 k 2 . The parameter s 2 controls to what 
degree the low energy modes are suppressed. To com- 
pletely exclude the effect of low energy modes, we should 
let s 2 — > oo in the end of calculations. 

Now we are to write down a typical regulator. An 
immediate choice may be a step-type regulator, due to 
its simplicity, 

r(q 2 /k 2 ; s 2 ) = (s 2 k 2 /q 2 - 1)6(1 - q 2 /k 2 ). (5) 

However this regulator is not continuous, and can raise 
problems. So we need for continuous ones. There is a lot 
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FIG. 2: Logistic-type regulators with different values of s 2 . 



of smooth approximations of the Heaviside step function 
Q(x), so we can use these functions to construct smooth 
regulators. For example, a logistic-type regulator, 



r(<z 2 A 2 ;s 2 ) 



i 



l 



\ _ e -(? 2 /(s 2 fe 2 ) ^ _|_ e -s 2 (l-g 2 /fe 2 ) ' 



(6) 



an erf- type regulator, 



1 + erf(s 2 ) 



and an arctan-type regulator, 



r(<r7*V) 



k 2 _ q i /k 2 tt/2 + arctan(s 2 (l - q 2 /k 2 )) 



tc/2 + arctan(s 2 



(8) 

In all cases, the second factor is a smooth approxi- 
mation of the step function 0(1 — q 2 /k 2 ) which is ar- 
rived at the limit s 2 — > oo. The first factor is to ensure 
r{q 2 /k 2 ]s 2 ) ~ s 2 k 2 /q 2 when q 2 < k 2 . When s 2 -t oo, 
for all above cases, r(q 2 /k 2 ; s 2 )-q 2 tend to an infinity step 
function which takes the value of infinity at q 2 < k 2 and 
the value of zero at q 2 > k 2 . One can see this trend in 
FIG. [21 where we have plotted profiles of the logistic- type 
regulator with different values of s 2 . 

Actually, the trend of tending to a step function and 
the trend of tending to infinite suppression in the left 
and no suppression in the right can always be parame- 
terized by two independent parameters. However, we use 
only one parameter s 2 , since our later discussions do not 
depend on this simplicity assumption. 

To implement numerical computations, we will take a 
folded-straightline regulator, in the form 



r( q 2 /k 2 ; S 2 )-q 2 /k 2 =^-^(^-l 



(9) 



in the interval [1 — e, 1 +e] , and constant elsewhere, s 2 — > 
oo is equivalent to £ — > oo and e — > 0. The order of 
taking limit is irrelevant. 



IV. APPLICATION TO QCD 

In this section we will apply the functional renormal- 
ization group method to QCD to study its renormaliza- 
tion flow. 

The Lagrangian of QCD in Euclidian spacetime is, [3l[ 
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(d"A a ) 2 + c a d"D ac c c 



+ c a d 2 c a + gf abc c a d ll (A b l c c ) 

+ gr bc (d K Al)A b K A c x + l -g 2 r he r dE A a K A\A^A d x . (10) 

We have used the convention {7^,7^} = —28^1^^ and 
the gauge £ = 1. 

Now go to momentum space and take the following 
ansatz for the effective action 



r = / Z^°'(p)(jJ + m)V> a/ (p) 
Jp 

+ ±Z A J AX-p)(p>8>»' - (1 - |K/)i:(p) 

- zj c a (p)p 2 c a (p) 
Jp 

+ gZ 4 ,^Z A \( [ rf(p)rt a A«(p-p')r f (p') 



+ 9 Z c Jzl / f abc (ip»)c a {p)A b ^p~p')c c {p') 
+ gZ A ^Z~ A [ [ [(2n) d S^(p + p' + q) 

J p J p' J q 

■ r bc (- iq ,)A a MA»Al{p') 



-a 2 7 2 

4 9 Z a 



(2ir) d 5( d \p + p' + q + q') 



'p Jp' Jq Jq' 

■ r be f^A a K {p)Ai{p')A c K (q)A d x (q'). (11) 

For simplicity, we have assumed all fermions have the 
same mass. 

The regulators are 

as,. : / / r f (p)R a / f (p^^ b, (^ ( 12 ) 



pdq 



1 



AN, = 2 JJ^(-p)R a A b ^(p,q)A b l/ (q); (13) 



AS C = / c a (p)R a c b (p,q)c b (q), 
Jp 



(14) 



with 



R a / bf (p,q) = Z^Spqr^/k 2 ;^ (15) 
R a A^(p,q) = Z A S ab S^S pq r A (p 2 /k 2 ;s 2 )p 2 , (16) 
Rf(p,q) = -Z c S ab S pq r A (p 2 /k 2 ; S 2 )p 2 7 (17) 
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with (1 + t>) 2 = (l + r A ). 

The functional renormalization group equation can be 
derived as 



d t T = ^Tr[(d t R A )G AA ]-Tr[(d t R^)G^} 
- Tr[(d t R c )G cc ], 



(18) 



where Gaa, G^ and G c e are connected Green functions, 



G AA = (A(q)A(p))-(A(q))(A(p)), 

= MffW(p))-Wg)>$(p)>, 
G CE = (c(q)c(p))-(c(q))(c(p)). 



(19) 
(20) 
(21) 



Since G ■ f < 2 > = 1, one has 



G AA 



=,(2) _ 

^(2) 
AA 



f(2) f (2)-lf(2) 
1 ,/iA 1 1 Ai/> 



^(2)f(2)-lp(2) _ p(2)f,(2)-lp(2) 
Ac cc 1 C A cc 1 cA 



f(2)f(2)-lf,(2) _ p(2)p(2)-lp(2) 

1 Aj/> tpijj ijjA Aip tfiip ipA 



G r ?i — 



=>(2) 



f(2)f(2)-lf,(2) 
cA AA 1 Ac 



(22) 

",(23) 
(24) 



After a lengthy calculation, we arrived at the following 
flow equations for the coupling g and the mass m (see 
Appendix [Al for the details) 



= ~g C 2 (r)K m , 

m 

dtg 



(25) 



-j = -g 2 [C 2 (r)K 1 +C 2 (G)K 2 -N f C(r)K 3 ](26) 

Dimensionless quantities are defined as m = m/k, g 2 = 
[Jdn d /(2Tr) d }g 2 /k 4 - d and 



K m = d[J(l,4,l;m 2 ) + J(0,2,2;m 2 )] 
(d- l)(d- 2) 



Here J(a, 6, c; m 2 ) are dimensionless momentum integrals 
defined by 



J(o, 6, c; m ) 



(«9 t rA) • fc 2 ( a+c )- d (27r) d //^<i 
, (/ 2 ) a (l + r^) & [m 2 + ; 2 (l + r^) 2 ]' 



.(27) 



This integral encodes all regulator-relating informa- 
tion. It will be investigated closely in the next section. 



V. MOMENTUM INTEGRAL AND 
REGULATOR-INDEPENDENCE 

Since the Wetterich equation has a one-loop structure, 
momentum integrals appearing in flow equations all have 
the form of (|27p , provided there is only one mass param- 
eter in the considered theory. 

Generally the following relations hold, 



J(a, b, c; m 2 ) 

= to 2 J(a, b, c + 1; m 2 ) + J(a — 1, b — 2, c + 1; m 



dJ(a, 6, c; m 2 ) 
dm 2 



-cJ(a, b, c+ 1; m 2 ). 



(29) 



First consider the integrals with a + c = d/2. The 
integration can be easily carried out analytically since 
kdkr A = —ldir A . When m 2 = and 6/2 + c - 1 > 0, the 
result is 



J(a, 6, c; 0) 



1 



c- 1 



(30) 



When m 2 G [0, 1) and 6/2 + c — 1 > 0, the result can be 
given by the series, 

CO , , 

J(a, 6, c; to 2 ) = V ^ J(a, 6, c + n; 0)(-to 2 )™, (31) 

n=0 



K 2 = - 



where (c)„ is the rising factorial, defined as (c) = 1, 
(c)i = c, (c)„ = c(c + l)—(c + n — 1). In this situation, 
since the momentum integral can be carried out analyt- 
ically without any reference to the concrete form of the 
4(d — 2) ^ g 3- to 2 ) — ~ l)(d ~ 2) g ^ ^regulator, these results are regulator-independent. 

d d > > > ! integrals m £h e fj ow equations (1231) and (|2l))) are 

all of this kind when d = 4. When to 2 = 0, the K 
coefficients can be calculated as: 



-J(1,3,1;to 2 ); 

Ki = (d- 2)[J(1,4,1;to 2 ) + J(0,2,2;to 2 )] 
2 



d-2 



[J(1,4,1;to 2 ) + J(0,2,2;m 2 )] 



+ ^4j(-l,0,3;m 2 ) 



d 

3(d- 1) 



if m = 3, = 0, X 2 = 



d 

16(d-2) 



J(l, 5, 1; to 2 ) + ^Lil j( , 3, 2; to 2 ) 



11 

3x2' 



3x2 



d(d + 2) 



d- 14 8 
2 + d 



J(2,6,0;to 2 ); 



^^-^(-l ) -2 ) 3 ; ^) + ^±| J (-2,-4,4;^ ) 



G4 



rf(d + 2) 



J(-3,-6,5;m 2 ^ 



The ordinary one-loop /3 function of QCD and the ordi- 
nary mass running equation are reproduced. 

The integrals with a + c d/2 are little more compli- 
cated, though generally the integration can not be carried 
out analytically. As we will eventually take limit after in- 
tegration, we can make an approximation for the integral. 
As has been aforementioned, the flow is driven only by 
modes in the vicinity of the cutoff scale, as indicated by 
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the fact that dtTA peaks at q 2 /k 2 = 1. Actually the peak 
of the integrand is rather sharp, and values at elsewhere 
are suppressed to nearly zero if s 2 is sufficiently large. So 
we can make the approximation: 



dl 2 ■ integrand 



l+e 



dl 2 ■ integrand, 



(32) 



hen s — > oo, e — > 0. 

The integral J(a, 6, c; 0) can be evaluated as: 



z 2 =i-< 



l d - x dl 



d t r A ■ fc 2 (°+ c )- d 

(p)a+c( 1+rj4 )6/2+ C 



(33) 



When a + c - d/2 = and 6/2 + c - 1 > 0, this integral 
can be carried out directly, yielding 1/(6/2 + c — 1) if 
,s 2 — > oo in the end, which is the same result as derived 
formerly. When a + c - d/2 ^ and 6/2 + c - 1 > 0, this 
integral can also be done through the mean value theorem 
of integration, yielding the same result 1/(6/2 + c — 1) if 
s — > oo in the end. These results are both independent 
of the specific form of the regulator rA(q 2 /k 2 ; s 2 ), hence 
regulator- independent. 



VI. FLOW PATTERNS AND CHIRAL 
SYMMETRY BREAKING 

In this section, we first explore different kinds of flow 
patterns. A mechanism for the restoration of chiral sym- 
metry will be arrived at. Then we prove this mechanism 
indeed applies to QCD. In this section the discussions 
will be constrained to situations of d = 4. 

First consider the flow patterns of the one-loop renor- 
malization group of QCD, (|25t and ([26)) . The coeffi- 
cient of g 2 in the equation (1261) always has a zero only if 
Nf > 11 x 3/2. While K m is always positive, so there is 
no non-trivial fixed-point for the flow equations (1251) and 
(1261) . Shown in FIG. [3] are their flow patterns. The left 
panel is typical for Nf < 11 x 3/2, and the right panel is 
typical for N f > 11 x 3/2. 

These flow patterns do not capture the non-occurrence 
of dynamical chiral symmetry breaking. However there 
are other kinds of flow patterns that may do so. To ex- 
plore these ones, it is natural to add in two loop contri- 
butions as the next step. 

Let's recall the ordinary two-loop renormalization 
group equation for the mass (32| . 




m'lk„ 



N=10 




FIG. 3: Flow patterns with d = 4 and complete suppression. 
The arrows point to infrared. The red dashed line stands for 
/3 = 0. 



flow equation of the mass will develop a zero point. In 
turn, the m 2 -g 2 flows will develop a nontrivial fixed point. 
However, coefficients in this equation are independent of 
m 2 . In other words, this equation corresponds to Eq. 
(|2~51) with to 2 = but at two-loop. If we use it directly 
to depict the flows, the fixed point will be displaced and 
streamlines distorted. Even worse, the flow patterns may 
be changed. Nevertheless, to see how flow patterns cap- 
ture the non-occurrence of dynamical chiral symmetry 
breaking, it will be enough to make a parametrization of 
Eq. 

Literally, we multiply two-loop terms in Eq. (|34[) by 
factors /(m 2 + f3i) ai (a* and fti are constant param- 
eters), and add them to Eq. (I2~5"j) . This form of mass 
dependence is motivated by the J- integral (|2T|) . ai is 
roughly the number of fermion propagators in two loop 
momentum integrals. The difference among fii can be 
generated after loop momentum integration. 

The resulting equation of these manipulations will be 
used to explore different flow patterns. One can see that 
the resulting equation may not be the true one of QCD, 
but it is enough for qualitative investigation. 

Concretely, the following system will be considered: 



dtfh 2 
2m 2 

2g 2 



l+g 2 C 2 {r)K m + H 2 (rh 2 ){g 2 ) 2 



(37) 



-g 2 [C 2 (r)X! + C 2 (G)K 2 - N f C(r)K 3 p8) 



with 



dtfh 2 



-2m 2 (l 



hig 2 



h 2 (?n 



(34) 



H 2 {m 2 ) = h 2 \p°^ + h 22 p2 2 - h 23 pz 



(39) 



with 



where 



hi 
h 2 



3C 2 (r), 

3 Q7 

7 C 2 (r) 2 + -C 2 (r)C 2 (G) 
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(35) 

C 2 (r)N f . (36) 



Observe that, in h 2 the coefficient before Nf is negative. 
So for large enough Nf , this term will dominate and the 
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h 2 i = -C 2 {r)C 2 {G), 
h 22 = \c 2 {r) 2 , 
h 23 = \c 2 {r)N f , 



(40) 
(41) 
(42) 



6 



Pi 

m 2 + Pi 



1,2,3. 



(43) 



This system is not the one of QCD. However, since 
our purpose is to explore different kinds of flow patterns 
rather than to discuss QCD directly, it is reasonable to 
use such a system. Actually, to explore different kinds of 
flow patterns, any functions are permissible in the right 
hand of the flow equations (|3"T|) and (I3"51) . The only prob- 
lem is whether they can lead to new and meaningful flow 
patterns. 

One loop flow equations of QCD, (1231) and (|2"o| . do 
not capture the restoration of chiral symmetry, but they 
provide a base to be perturbed and to generate new flow 
equations. 

Besides the trivial fixed point (m 2 = 0,<? 2 = 0), the 
flow equations (|3T[) and (|38[) now may have a nontrivial 
fixed point. Taking m 2 as the argument, the possible 
nontrivial fixed point is 



at; = 



C 2 {r)K 1 + C 2 {G)K 2 

C(r)K 3 
-C 2 {r)K m - y/[C 2 (r)K n 
2H 2 



(44) 
(45) 



Another unphysical root has been discarded. To ensure 
g 2 * > 0, the sufficient and necessary condition is: 



H 2 (m 2 *) < 0. 



(46) 



To see the behavior of flows near the above fixed point, 
we estimate the critical exponents. If Ai and A2 are the 
critical exponents of this fixed point, then, 



A1A2 



Ai + A 2 



-1/g 2 * + H 2 (m 2 



>0; 



-2m 2 * 

-2to 2 * 
dH 2 



~g 2 *C 2 (r) 
dH 2 



dK n 



din 2 



dm 2 

(g 2 r 



+ 



dHi 



dm 2 



(9 



~2*\2 



dm 2 



When 



dHo 



dm 2 



<0 (>0), 



(47) 



the fixed point is IR-attractive (IR- repulsive). 

One sees that the fixed-point structure of Eqs. 
and (|3"5)) depends on the choice of parameters (on, Pi). 
For example, the choice of (a>i, 012,0.3, fix, P2>Ps) — 
(3,4,5,1,1,36) will lead to an IR-attractive nontrivial 
fixed point for all Nf > 17.12. While the choice of 
(2,3,4,1,2,5) will result in an IR-repulsive fixed point 
for all Nf > 25.56. We have depicted some flow patterns 
in FIG. S|- FIG. El 




FIG. 4: IR-attractive nontrivial fixed point. The arrows point 
to infrared, (ai, Pi) = (3, 4, 5, 1, 1, 36). 




FIG. 5: No nontrivial fixed point. The arrows point to in- 
frared. {a.i,Pi) = (3,4,5,1,1,36). 




srW) 



FIG. 6: IR-replusive nontrivial fixed point. The arrows point 
to infrared, (on, Pi) = (2, 3, 4, 1, 2, 5). 



7 



Shown in FIG. [4] is the flow pattern for Nf — 18 with 
(ai,/3i) — (3,4,5,1,1,36). One sees that all flows are 
going to the IR-attractive fixed point. At this point, 
to 2 is a constant, which means to 2 = m 2 k 2 goes to zero 
when k 2 — > 0. So there is no mass generation or chiral 
symmetry breaking. 

Dynamical chiral symmetry breaking is in principle a 
non-perurbative behavior, much the same as color con- 
finement. However, color confinement can be gauged 
via perturbative calculation, in the form ol asymptotic 
freedom. Similarly, dynamical chiral symmetry breaking 
may also be gauged via perturbative calculation. To chi- 
ral symmetry the fermion mass is more relevant than the 
gauge coupling. So the effect of fermion mass seems more 
important than higher order quantum corrections. Thus 
the above discussion from the perspective of perturbative 
renormaliztion flow is meaningful. 

Shown in FIG. [5] is the flow pattern for Nf — 10 with 
(aj,/3i) = (3,4,5,1,1,36). There is no nontrivial fixed 
point and all flows to the limit (oo, m 2 (H 2 = 0))- This 
cannot be interpreted as a restoration of chiral symmetry, 
since the coupling g 2 has gone beyond the perturbative 
region. Shown in FIG. [6] is the flow pattern for Nf — 30 
with (ctfj, /3j) = (2, 3, 4, 1, 2, 5). An IR-repulsive nontrivial 
fixed point turns up. 

In the above, different flow patterns are displayed. One 
of them, in FIG. 21 provides a mechanism to explain the 
non-occurrence of dynamical chiral symmetry breaking. 
Theories attracted to an IR-attractive nontrivial fixed 
point of a finite to 2 * will show no dynamical chiral sym- 
metry breaking. 

Now we show this mechanism can also apply to QCD. 
We shall prove the inverse proposition: a non-occurrence 
of dynamical chiral symmetry breaking indicates the the- 
ory is attracted to an IR-attractive nontrivial fixed point 
of a finite to 2 *. 

The mass running equation of all-loop QCD can be 
written as 



dtfh 2 



-2m 2 



Here Hi(m 2 ; Nf) is the mass-dependent coefficient of the 
i-loop correction to the fermion mass in QCD. Restora- 
tion of chiral symmetry requires m 2 = in the infrared. 
So the scaling exponent of to 2 is smaller than 2, or the 
value of sum of terms in the square bracket is smaller than 
1. This is possible, because when Nf large enough, some 
of Hi(m 2 ;Nf) will be negative. Since Hi(m 2 ;Nf) are 
loop corrections to the fermion mass, they will all tend 
to when m 2 —> oo because of the fermion propagators 
in the loop integrals. Actually when m 2 large enough 
and g 2 small enough, the sum of quantum corrections 
will be positive. This is due to the fact: loop integrals 
with more than one fermion propagators then will sub- 
dominate, and the one-loop contribution will dominate. 

This discussion can also apply to the non-perturbative 
situation. Loop integrals with more than one fermion 



propagators sub-dominate only if to 2 is sufficiently large. 
The difference is that the one-loop integral is now with 
a gluon propagator that has been renormalized in pure 
gauge theory. It is still positive, for a gluon cannot 
change to a ghost by mere renormalization. One can see 
that the above discussion is not based on perturbative 
expansion, but on 1/to 2 expansion. 

So, if there is no dynamical chiral symmetry breaking, 
to 2 should be a finite value in the infrared. The coupling 
g 2 should be a nonzero (finite) value as well. The theory 
is attracted to an IR-attractive nontrivial fixed point of 
a finite to 2 *. Thus we proved the proposition. 



VII. THE CONFORMAL WINDOW 

In this section, we will see how conformal window is 
reflected in renormalization flows. 

Recall the ordinary two-loop running equation of the 
coupling in QCD: 



with 



11 



-C 2 (G) - -C(r)N f , 



(48) 



(49) 



3 

fa = l ^C 2 {G) 2 - l ^C 2 {G)C{r)N } -2C 2 (r)C{r)N^) 

One can do similar manipulations as those for the mass 
running equation, however there is another choice which 
is simpler and more certain. 

In Eq. (|48p . there are two turning points for Nf at 
which the behavior of (3p(= dtg 2 ) changes. (3 g 2 is always 
positive for Nf > 11 x 3/2 and always negative for Nf < 
51 x 3/19. When 51 x 3/19 < N f < 11 x 3/2, the sign of 
changes at a nontrivial fixed point, g 2 * = — fi/ fa. 

This property can also be found in the one-loop 0- 
function of d-dimensional QCD with d = 4 — e < 4: 



d t g 2 = -g 2 (e + g 2 F 1 (0)), 



(51) 



where iq(0) is the value of 2[C 2 (r)K 1 + C 2 {G)K 2 - 
N f C(r)K 3 ] with to 2 = 0. When N f < \C 2 (r)K x + 
C 2 (G)K 2 }/[C(r)K 3 ] at to 2 = 0, /3 g 2 is always neg- 
ative. When Nf > [C 2 (r)Ki + C 2 {G)K 2 ]/[C{r)K 3 ] 
at to 2 = 0, the sign of j3p changes at a nontrivial 
fixed point, g 2 * = —e/Fi(0). The difference is that 
[C 2 (r)Ki + C 2 (G)K 2 ]/[C(r)K 3 ] now plays the role of the 
lower turning point, while the larger turning point now 
goes to infinity. 

One can see this similarity from FIG. [7] So we can use 
the one-loop /3-function of d-dimensional QCD to mimic 
the two- loop /3-function of 4-dimensional QCD. The fol- 
lowing system will be considered: 



<9 t TO 2 = -2to 2 



(52) 



l + ~g 2 C 2 (r)K m + H 2 (rh 2 ){g 2 ) 2 

d t g 2 = -2(g 2 ) 2 [C 2 (r)K 1 +C 2 (G)K 2 -NfC(r)K 3 ] 

(53) 
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+ 


11x3/2 


[ + 






51X3/19 




Two-loop QCD in 4 dim 



2 4 6 8 10 



One-loop QCD in 3.5 dim 



1 2 3 4 5 6 7 

\tilde{g}"2 



FIG. 7: The similarity between the two situations. 




FIG. 8: IR-attractive nontrivial fixed points. The arrows 
point to infrared, (on, pi) = (3, 4, 5, 1, 1, 36). 



This system deviates slightly from the one of the last 
section, so the flows will not be changed drastically. Since 
now there is a nontrivial fixed point on the axis of the 
coupling, there would be modifications. 

Shown in FIG. [8] is the flow pattern for the situation 
of d = 3.5, N f = 24 and (o^A) = (3,4,5,1,1,36). The 
pattern is much the same as that of FIG. 2] except for 
the fixed point on the axis. All flows are going to the 
nontrivial fixed point with a finite dimensionless mass, 
which means there is still no chiral symmetry breaking. 
Mapped onto two-loop d = 4 QCD, there is no chiral 
symmetry breaking when Nf < 11x3/2 and larger than 
some critical number. So a conformal window opens up. 



group method to theories with both fermions and bosons. 
Hopefully our calculation can serve to help clarify. 

A parameter was introduced to discuss the dependence 
of flow equations on regulators and it was shown that 
flow equations will be independent of regulators if this 
parameter is sent to infinity in the end of calculation. 

Different flow patterns in the m 2 -g 2 plane were ex- 
plored. (By no means, all possible ones were reached.) 
One of them provides a mechanism to explain non- 
occurrence of dynamical chiral symmetry breaking. How- 
ever this flow pattern was obtained from flow equations 
that were not of QCD. To show the relevance of this 
mechanism to QCD, it was proved that it indeed applies 
to QCD, in terms of an 1/m 2 expansion of the fermion 
mass renormalization group equation of QCD. One saw 
that dynamical chiral symmetry breaking cannot occur 
if and only if the theory is attracted to an IR-attractive 
non-trivial fixed point of a finite dimensionless mass. 

The existence of the conformal window was also dis- 
cussed in the language of renormalization flow. 
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Appendix A: Details 



1. Preparation 



If we label the fields as (ip, ip, A, c, c), and let 
(Jj,, Jjf,, Ja, Jc, Jc) be their sources, and W be the gener- 
ation functional of connected Green functions, the two- 
point connected Green function has the structure: 



G 



C - 
^tp^p 



\ 



G 



A4< 



\ 



G^A 

Ga%P GaA G Ac Gac 

GcA ■■■ Gcc 

G C A G C c 



(Al) 



The terms represented by ellipsis will not be used. Other 
terms are defined by 



S 2 W 



G - = 



G 



TpA 



S 2 W 
5 L Jj,SJ A 



(A2) 



The r Hessian matrix has the structure, 



VIII. CONCLUSION 

In conclusion, the functional renormalization group 
method was applied to QCD. There are variations in lit- 
erature about application of functional renormalization 



r (2) = 



/ o 

r (2) 

M 

Aip 





r (2) 



r (2) 
1 Aip 






F (2) 

ip A 
,(2) 
ipA 
.(2) 
AA 
.(2) 
cA 
.(2) 
iA 



r ( 



o 
o 

1 Ac 



r( 2 ) 






- Ac 
.(2) 

- cc 

o J 



(A3) 



Elements of the T Hessian matrix are, 
S 2 T 



r (2) = 



,(2) 



r (2) _ 



6 L c a (l)6 R c b (l') 
= -Z c 5 ab 5 w l 2 - gZ c ^r bc (U K )Al(l - I 1 )- 

8 2 T 



S L c b (l')S R c-(l) 

Z c 5 ab S w l 2 + gZ c ^Z^r bc (U K )A c K (l - I')- 
5 2 T 



r (2) 

L Ac 



' ' ~ 5 L c-(l)5A b (l>) 
5 2 T 



5Al{-l')8 R c a {l) 
S 2 T 



-=gZ c ^Z^r bc (a u )c c (l-l'); 

-gZ c ^Z A ~t bc {il v )c c {l + 0; 



F (2) _ 
L Ac ~ 



6A-(-l)5 R c»(l>) 

?cab- 



.(2) 
- cA 



gZ c ^/Z A f cab i(-l^ + l'^)c c (-l + I')- 
6 2 T 

' d L c b (l')SA a U) 



--(2) 



= -gZ c ^Z A f cab i(l ll + 1^(1 + 1'); 



5 2 r 



W 5 L ip a f(l)S R ^ b f(l') 

= Z i ,5 a ' b f6 w (jl + m) + gZ^^Z~ A b a s b s 1 n a Al(l - I'); 

r ( 2) s 2 r 

^ 6 L ^ bf (l')5 R ip a f(l) 

= -Z^S a f b fS ll ,f + m)-gZ^^Z^S^ b nrt a ] T A-(l 

where T means transposition; 

6 2 T 



+ gZ A ^Z A r hc ^\-i){l K + QA c K (l - V) 
+ g 2 Z 2 A [ A^{q)A dv {l-l' - q) 

Jq 

^ jeab jecd jead j e ^ c ^ 

+ g 2 Z\ f A cX (q)A dX (l-l'-q)S^ 

Jq 

}_^cac jcbd _|_ jead jebc^ 

Derivative expansions of the regulated propagators are 
1 



1 

~ {Jt+i>){l + ri,{l 2 lk 2 )) + m 
_ m-/(l + r^,) + 



to 2 + P(1 + t>) 2 m 2 + P(l + r^,) 2 
-2p-i(l + r v ,) 2 



+ [m-/(l + r>)] 
+ [m-/(l + ?>)] 



[m 2 + / 2 (l + r^) 2 ] 2 

-p 2 (l+r>) 2 
[m 2 + / 2 (l + r^,) 2 ] 2 
2p-Z(l + 7>) 2 



212 



[m 2 + Z 2 (l + r>) 



-0, 



(l+p) 2 (l + r A ((l+p) 2 /k 2 )) 
1 

~ G+rt 2 (l + ^(/ 2 /fc 2 )) 
1 2p- 1 



l 2 (l + r A ) (l 2 ) 2 (l + r A ) 



,(2) _ 

* A <5i,V a ^(/)^(/') 

£ 2 r 



.(2) 



,(2) 



„± l? , + {p % (A5) 



(Z 2 ) 2 (l + r A ) (l 2 ) 3 (l + r A ) 



The d t derivatives of the regulated propagators are 



8 2 T 



r ( 2) = s 2 r 

+ A W(i')<*4£(0 
(J 2 ! 1 



-=gZ^ v fz A i> b f(-l + l')rt a ; 
gZ^^Z A ~[i, b f{l + l') 1 H a ] T ] 



/(l + r,/,)+m v /(l + r^,)+m 

[m 2 -Z 2 (l + rv) 2 y + 2mZ 2 (l + r>) 



= (<9 t r^)- 
= Fi/ + F 2 m. 



+ Z 2 (l + r^) 2 ] 2 



.(2) 
AA 



Z A 5 ab 6^6 u ,l 2 



SA-(-l)6A b (l') 



+ gZ A \J Z A f i{2l u — l v )A c ^(l - I') 
+ gZ A ^Z~ A r bc i- % ){l ll -2l' tl )Al(l-l') 



' -{d t r A )l 2 ' 



l 2 (l + r A ) y l 2 (l + r A ) 



EE B. 



(A6) 



(A7) 
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2. Fermion propagator 



To get the flow equation for the fermion propagator, we encounter 

1 



g'Z^ZA^i-l + n^t 11 ^ 



/"(l + r v ,)+m' ' r v " w Z A /' 2 (l + r A ) 



/" T (l + r^) + m l 



/,{! + r^)+m 



g 2 z 2 z Al n a r'{i-i" 



g 2 Z^C 2 (r)(-d) / ^°/(p)^/(p)m. 



Z A l"\l + r A ) 
B F 2 



ip b f(-l" + l')^t a -j 



1 Z A l' 2 (l + r A ) 
l/Z^ 



/ (1 + r>) + m 



+ ff %C- 2 (r)(d-2) / r'(p)W(p) 
+ higher order terms. 
The flow equations for m and Z^ can be read off as 



TO 2 +Z 2 (l + r>) 2 Z 2 (l + r A ). 

B(1 + t>) 2B(l + »v) 3 i 2 /d 2Fi/d 1 



+ Z 2 (l + r^) 2 + [m 2 + Z 2 (l + r^,) 2 ] 2 + / 2 (1 + r A ) J 



5 t™ ~2r< i \ 

= -g C 2 (r) 

to 

dtZ^p „ 2 i \ 

= -g C 2 (r) 



dJ(l, 4, 1; to 2 ) + dJ{0, 2, 2; m 2 ) - ^ ^ 2) J(1, 3, 1; to 2 ) 



(d- l)(d-2) 

d 



J(l,3,l;m 2 ) 



(A8) 

(A9) 
(A10) 



3. Gluon propagator 



Here one encounters 



Tr(2U&r A Z 2 ) 



Z A Z 2 (l+r A ) 



9 2 Z\{...} 



Z A l' 2 (l + r A ) 



0. 



Here g 2 Z A {...} indicates terms in T^ A that come from the 4-gluon vertex. 
To get the flow equation for the gluon propagator, we encounter 



Z A l"H\ + —f Z ^i-} zj^l + r A ) 



-Tr(- Z c d t r A l 2 ) 



-TriZ^dtr^) 



-1 



z c p(i + r A y 



gZ c ^/Z A r e {il K )Al{l^l") 



-1 



l/Z^ 



/(l + ?>) + to" 



Z c /" 2 (l + r A ) 



gZcy/zIf^'^AiiV'-l') 



-1 



/"(l + r>)+' 



Z c /' 2 (l + r A ) 
1/^ 



/ (1 + ?>) +m 



= g 2 Z A C 2 {G) J^{- P ){5^p 2 -p^ v )Al{p)- J i 
-g 2 Z A N f C(r) J^A;(-p)(S^p 2 - Pfl p v )Al(p) ■ jf{. 



« 2 (l + r A ) 



16d-32 d-14 



d(d + 2) 2 

^(1 + ry,) 3 ? 2 -16 F!(l+^) 5 g 2 ) 2 64 ) 
[TO 2 + / 2 (l + r^) 2 ] 2 d + [TO 2 + / 2 (l + r^) 2 ]3d(d + 2)J 
+ higher order terms and gauge symmetry breaking terms. (All) 



Here gZ A ^/Z A {...} indicates terms in T^ A that come from the 3-gluon vertex. 
The flow equation for Z A can be read off as 



l d t Z A 
2 Z A 



= -~9 2 C 2 {G) 



- ~g A N f C{r) 



16d-32 d-14 8 
d(d + 2) + 2 + d 



J(2,6,0;to 2 ) 



64 



J(— 3, —6, 5; fa 



. (A12) 
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4. Fermion-gluon vertex 

The flow equation for the coupling can be derived from the renormalized fermion-gluon vertex. 
1 -MZ A d t r A l 2 )——^- i -g 2 Z^Z A ^(-l + /"h^ 77777 ffV^W - l"')Yt e 



2 'Z A P(l + r A y w ^ v " /"(l + r^+rn " py " /"'(l+r^+m 



Z A /' 2 (l + r^) 



■[f'(i / + i'")/t a,T ! 



Z A Z' 2 (l + r A ) 



'/(l + r^J+m" ^ - P v /"(l + ^) +m 

-Tr(Z^ t r^/) /(1 + ^ + m (- 9 Z^y/z2)^r f Q - O + ~ } ( 9 Z^y/z2)^ (-1" + l'")-ft a 



±¥± gZtyfzlA'Al'" - l')Yt e -r 



1/Zj, 



/'"(l + r>)+rn " v /'(l + rv,)+m 

+ ^Tr(Z A 9 t r^ 2 ) ) {-g 2 Z^^)^ {-I + l")^ 1 7^"' (*" ~ H 

2 ZA<^(l+rA) / (l + r^,)+m 

- — z7{. ..};■-/■ l — 



Z A l"»(l + r A y " v Z A Z'2(l + r A ) 

+ l^ r ^ Z A P(l\rx) gZ ^ { --->^ zW + rx) 

• (-s%V^W 'H" + i w )7 p *° 77777: 1 : 7 ^V a/ ('"' - 0" 



/"(l + r v ,)+m' " v ^Z' 2 (1 + t-a) 

-i-R(Z A 9 t r A f 2 ) 1 (- ff 2 Z v , A /zI)[tV^a ,/ + /)] r ,^, 1 \ $*> (l'" + l")-ftf 

A Zj A 1 (L+r A ) j[ (1 + r,p) + m 



Z A l>"2(l + r A y " v ^-""'i' Z A l>2(l + r A ) 

1 m t w-r c\ i9\ 1 r , I r-, r -1 nr 11 n 1 



• (-^^)[tYf / r+i' , f -^- 1 - — [^a' + r) 7 v t 6 ] r - ! 



/"^(l + r^+m 1 " V ; ' J ^' 2 (l+r A ) 



/(l + r^) + m v * vv r v ; ^" 2 (l + r A ) a vv ^ L - J << Z A V^{l + r A ) 



■gZ^y/z2^{J,"' -/')7^ a 777 



1/^ 



/ (1 + »>) +m 
g 3 Z^^ [ [ $ a '(p)'ffAl(p-p')1> a, (p') 

J p J p' 

l 2 (l + r^) 2 {d - 2) 2 /d + m 2 {d - 2) 



- [C 2 (r)-lc 2 (G)]^B 
2[C 2 (r)-ic 2 (G)]^ 



[to 2 + P(1 + t>) 2 ] 2 
-Fi? 2 (l + ry)(rf-2) 2 /rf + F 2 m 2 (d-2) 1 

m 2 + F(l + rv,) 2 l 2 (l + r A ) 



r(r] f B 4Z 2 (1 + T> ) d-l l r(n f 1 f 4/ 2 F /-l 

" C2(G) / P(l + r A )m 2 + / 2 (l + r„) 2 d ~ 2^ (G) / [P(l + r A )Y ^ 
higher order terms. (A13) 
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Then the flow equation for gZ^\J Za can be read off as 



(d - 2) J(l, 4, 1; irt j + {d- 2) J(0, 2, 2; m 2 ) 



^^J(-l,0,3;m 2 ) 



^(G) 



— — J(l,5,l;m 2 ) + — — J(0, 3, 2; m 2 ) 



(A14) 



Appendix B: Some identities 

With the convention {7^, 7"} = -2S^I 4x4 and d- 
dimensional Euclidian spacetimc, 



7V = (-d)W 



yxyy* = ( d _ 2 ) 7 p. 



// — (— 1 2 )Iaxa- 



ji-ypjt = z V - yip. 



(Bl) 



(B2) 



777 P /«7 M = (d - 2)/ 2 7 p - (d - 2)2/Z p . (B3) 



(B4) 



(B5) 



Tr( 7 'V) = -45" v . 
Color matrix products: 

t a t b t a = [C 2 (r) - ic 2 (G)]t b 



1 



^yrasc = -iC 2 (G)t c . 

The following expansion will be used frequently, 
1 



1 11 111 

~A + ~e ~A~ A 6 A + A € A € A 



(B6) 



(B7) 



(B8) 



(B9) 
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